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ABSTRACT 


The  problem  of  estimating  frequencies  of  sinusoids  buried  in  noise  has  been  of 
great  interest  in  both  military  and  civilian  applications.  In  particular,  in  Control  Systems 
with  flexible  appendages  the  sinusoidal  vibrations  can  cause  instabilities  and  degrade  the 
performance  of  the  overall  system.  In  this  thesis  the  problem  of  identifying  frequencies  of 
disturbances  in  flexible  systems  using  advanced  Digital  Signal  Processing  techniques 
such  as  filter  banks  and  Quadrature  Mirror  Filters  is  addressed.  In  a  number  of  situations 
there  is  a  need  to  design  a  controller  for  a  system  with  flexible  modes.  In  particular,  in 
space  applications  solar  panels  and  robotic  arms  introduce  flexible  modes  in  the  system 
which  degrades  the  performance.  In  these  kinds  of  applications,  the  frequencies  of  the 
flexible  modes  can  not  be  modeled  accurately  a  priori  and  they  can  change  according  to 
the  operating  conditions.  The  proposed  approach  is  tested  by  computer  simulations. 
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EXECUTIVE  SUMMARY 


Satellite  systems  are  being  deployed  on  a  regular  basis  these  days  with  various 
missions,  namely  civilian  or  military.  Among  them,  at  a  great  percentage,  are 
communications  and  observation.  However,  it  has  been  observed  that  for  these  cases  the 
existing  vibrations  from  the  flexible  surfaces  of  the  satellite  and  the  atmospheric 
interferences  add  colored  noise  and  consequently  undennine  the  system’s  perfonnance.  In 
this  thesis  an  approach  to  estimate  frequencies  in  noise  based  on  Filter  Banks  and 
Quadrature  Mirror  Filters  is  developed.  The  general  idea  is  to  decompose  the  signal  into 
components  at  various  frequency  bands  and  compute  the  relative  energy  in  the  single 
bands. 

The  Quadrature  Mirror  Filters  approach  is  based  on  the  Daubechies  filters  which 
guarantee  orthogonality  of  the  frequency  components.  The  advantage  of  this  approach  is 
that  it  can  be  suitable  for  a  non-uniform  decomposition  of  the  frequency  axis,  so  that  one 
can  have  a  higher  frequency  resolution  at  desired  parts  of  the  spectrum. 

The  Discrete  Fourier  Transfonn  (DFT)  Filter  Bank  approach  is  based  on  an 
efficient  bank  of  passband  filters  all  derived  by  modulation  of  a  prototype  lowpass  filter. 
In  this  case,  orthogonality  of  the  various  frequency  components  can  be  guaranteed  only 
by  the  use  of  ideal  filters,  which  clearly  cannot  be  achieved  in  practice.  This  can  be 
approximated  using  a  Finite  Impulse  Response  (FIR)  prototype  lowpass  filter  with 
sufficient  attenuation  in  the  stopband.  The  technique  is  applied  to  the  particular  case  of 
estimating  the  underdamped  frequencies  in  flexible  structures.  For  the  Quadrature  Mirror 
Filter  case  a  recursive  structure  is  devised,  which  provides  recursive  decomposition  after 
frequency  bands  with  higher  energy. 
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I.  BACKGROUND  THEORY  /  LITERATURE  REVIEW 


This  chapter  covers  the  theory  of  the  DFT  Filter  Banks  and  the  maximally 
decimated  DFT  Filter  Banks  used  in  this  thesis.  Furthermore,  the  Quadrature  Mirror 
Filters,  Daubechies  filters,  and  ideal  FIR  Filters,  which  constitute  the  theoretical  basis  for 
this  research,  are  covered  as  well. 

A.  DFT  FILTER  BANKS  /  TRANSMULTIPLEXERS 
1.  Introduction 

In  several  applications,  in  order  to  be  able  to  estimate  the  frequencies  of  sinusoids 
hidden  in  noise,  there  is  a  need  to  decompose  the  original  signal  into  its  subcomponents. 
This  process,  called  “analysis,”  is  the  first  part  in  the  signal  decomposition  with  “filter 
banks.”  The  other  part,  called  “synthesis,”  is  the  assembling  of  the  decomposed  signal 
again  into  one  signal.  The  exactly  reversed  procedure  is  performed  by  the 
“transmultiplexers”  and  will  be  addressed  later  on.  Based  on  [Cristil],  the  diagrams  in 
Figures  1  and  2  show  the  general  case  for  both  filter  banks  and  transmultiplexers. 


Figure  1.  Filter  Bank  [After:  Cristil] 


1 


Synthesis 


Analysis 


Figure  2.  Transmultiplexer  [After:  Cristil] 


The  benefit  of  using  Filter  Banks  for  frequency  estimation  is  twofold.  First,  with 
the  frequency  decomposition  it  can  be  determined  which  of  the  narrower  frequency  bands 
contains  the  higher  energy.  Second,  the  data  used  for  frequency  estimation  has  better 
noise  characteristics  in  terms  of  signal-to-noise  ratio  (SNR)  and  the  noise  is  more 
“white.” 


2.  DFT  Filter  Banks 

The  latter  property  is  particularly  attractive  when  one  wants  to  estimate  the 
frequency  using  a  model-based  approach.  For  the  design  of  the  filter  bank,  a  prototype 
lowpass  filter  H(co)  with  bandwidth  n !  M  is  used  and  from  that,  all  the  other  filters  are 
derived  by  shifting  the  filter  in  the  frequency  domain  as  given  by 

Hk(co)  =  H(co-A  (1.1) 

M 

for  k=0,...,M-l  as  shown  in  Figure  3.  By  doing  so,  no  spectrum  areas  are  left  uncovered 
and  the  bandwidth  is  not  exceeded. 
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Figure  3.  Frequency  Spectrum  Coverage  [After:  Cristil] 


It  is  easy  to  see  that  the  transfer  function  of  these  filters  can  be  written  as 
Hk(z)  =  H(e~Kk2*IM)z ) ,  k=0, 

If  H(z)  is  decomposed  into  its  noncausal  polyphase  components  as 


M- 1 


H(z)  =  YjzIE_,(zm),  l=0,...,M-l 

1= 0 

Equation  1.1  becomes 

M- 1 

Hk(z)  =  H(e-Kk2*,M)z )  =  ^  (e-J(k2«IM^\l 

M- 1 

_  ^  zl e~j(kl27r/M) g  e~j{k2xM / M ) 

1  for  all  integer  k. 


1=0 


’z)  E_,  ((e 

M- 1 


j{kl7T  I M)  \M 


z)  )  = 


)=Yjz'e 


j{kl2n / M)  t-t  /  M 


since  e 


/=o 

- j{k27tM I M)  • 


/=0 


If  the  aforementioned  equation  is  applied  for  k=0,  then 

H0(z)  =  E0(zM)  +  zE_1(zM)  +  z2E_^zM)  +  ...  +  zM-iEM_l(zM) 
which  can  be  written  in  matrix  fonn  as 


Hk(z)  =  [  1  1  ...  1] 


E0(z  ) 
zE_{(zm) 


zm-'E_m+1(zm) 


Similarly,  the  general  form  of  Equation  1.4  can  be  written  in  a  matrix  form  as 


(1.2) 

(1.3) 


(1.4) 


(1.5) 


(1.6) 
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- j{k2n!  M ) 


(1.7) 


HAZ)  =  \}  e 


E,{zM) 

zE_x(zm) 

zm~1E_m+1(zm) 


Combining  all  M  filters  in  one  vector,  the  matrix  becomes 


"  H0(z)  1 

"i 

i 

1 

1 

N 

1 _ 

Hx(z) 

= 

i 

-j(k2xlM) 

-j{k2x(M-\)IM) 

zE_x(zm) 

1 

a? 

i _ 

i 

g-  j  {k27t(M -1)  /  M) 

g-j(k2x(M-l)2 IM) 

_zm-1E_m+x(zm)_ 

(1.8) 


By  doing  this,  the  analysis  filter  bank  in  Figure  1  is  significantly  simplified  (less 
coefficients  per  filter)  and  can  be  implemented  as  in  Figure  4,  where  E^  are  the  polyphase 
components  of  the  prototype  filter. 


Figure  4.  Filter  Bank  Analysis  [From:  Cristil] 


In  this  thesis,  for  the  case  of  casual  FIR  filters,  the  fonnula  that  will  be  used  for 
the  polyphase  decomposition  is 

M- 1 

H(z)  =  YJz~'El(zM),  1=0,  (1.9) 

1=0 


4 


Due  to  the  available  block  for  unit  delay  in  Simulink  (Y7),  and  by  following  the  same 
procedure  as  before,  the  author  gets  the  matrix  below  and  the  implementation  of  Figure  5. 


"  Ho(Z)  1 
H,{z) 

'i  i 

1  ej(k2jr/M) 

1 

gj{k2jt(M -1)/ M) 

E0{zM) 

1 

'3 

_ i 

1  ej(k2n(M-l)!M) 

gj(k2x(M-\f/M) 

_z-M+lEM_x(zM)_ 

x[n] 


vM-i[n] 


vi[n] 


v0[n] 


Figure  5.  Thesis  Filter  Bank  Implementation 


3.  Maximally  Decimated  DFT  Filter  Banks 

The  analysis  filter  bank  is  based  on  a  number  (bank)  of  bandpass  filters  with 
shifted  frequency  bands  equally  arranged  to  cover  the  whole  spectrum.  The  maximally 
decimated  case  could  be  summarized  as  the  decomposition  of  a  signal  x[n]  into  M 
channels  Vo[n],  vi[n],...,  Vm-i[h],  with  each  channel  downsampled  by  M.  Similarly,  in  the 
reverse  process  of  synthesis,  each  of  the  M  signals  Vo[n],  vi[n],...,  Vm-i[ii]  are  upsampled 
before  being  combined  into  the  single  signal  y[n].  Figure  6  illustrates  a  network  with  a 
maximally  decimated  filter  bank  that  is  equivalent  with  the  network  in  Figure  2,  where 
L=M. 
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Analysis 


Synthesis 


Figure  6.  Maximally  Decimated  Filter  Bank  [From:  Cristil] 

Based  on  the  Noble  Identity  the  above  network  can  be  reconstructed  with  the 
downsamplers  repositioned  and  the  zM  of  the  filters  replaced  as  z.  Figure  7  shows  the 
implementation  of  the  analysis  and  synthesis  networks.  Since  the  two  operations  of  the 
DFT  and  the  IDFT  cancel  each  other  out,  the  outcome  is  the  equivalent  system  shown  in 
Figure  8.  Under  certain  conditions  on  the  prototype  filter,  perfect  reconstruction  can  be 
guaranteed. 
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For  the  case  of  perfect  reconstruction,  the  prototype  filter  H(z)  used  for  the  DFT 
analysis  network  and  the  prototype  filter  G(z)  used  for  the  DFT  synthesis  network  must 
be  such  that 

E_k(z)Fk(z)  =  ±-,k=0,...,M-l  (1-11) 


From  this  condition  it  can  be  seen  that  if  FIR  filters  are  used,  E.k(z)  and  Fk(z)  are  simply 
zero  degree  polynomials  (constant)  and  therefore  considerably  restrict  the  class  of  filters 
that  can  be  used.  In  fact,  it  is  easy  to  see  that  the  only  filters  that  can  be  used  in  the 
maximally  decimated  DFT  filter  banks  are  either  ideal  filters  or  FIR  filters  with  length  M. 
As  a  consequence,  to  guarantee  perfect  reconstruction,  i.e.,  no  loss  of  information,  there 
is  a  need  to  downsample  by  a  factor  smaller  than  M.  [Cristi],  [Cristil] 


B.  FIR  FILTERS 


The  Ideal  LP  filter  passes  frequencies  up  to  a  cutoff  frequency  coc  and  rejects  all 
frequencies  above  coc.  In  other  words,  its  frequency  response  is 

1,  0  <  co  <  0)„ 


HLP(eJ(°)  = 


0,  (o  <co<n 


(1.12) 


Since  its  impulse  response  has  an  infinite  length  and  it  is  non-causal  this  is  non- 
realizable  in  practice.  An  approximation  by  an  FIR  filter  can  be  parameterized  by  a 
passband  ripple  and  stopband  attenuation  as  shown  in  Figure  9. 
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Zero-phase  Response 


Figure  9.  Example  of  Approximation  of  an  Ideal  FIR  Filter  [From:  Losada] 


Finally,  it  can  be  shown  that  the  larger  the  order  of  the  filter,  the  smaller  the 
transition  width  and  the  ripples  will  be.  So,  depending  on  the  application  the  filters  are 
wanted  for,  the  decision  has  to  be  made  as  to  how  “ideal”  one  wants  them  to  be  by 
selecting  the  allowed  values  for  the  transition  region  as  well  as  the  passband  and 
stopband  ripples. 

C.  QUADRATURE  MIRROR  FILTERS 

Quadrature  Mirror  Filters  are  extensively  used  in  subband  coding  of  speech 
signal,  image  compression,  communication  systems,  etc.  With  these  filters,  an  input 
signal  can  be  decomposed  into  its  subband  components  and  then  resynthesized  by  them, 
with  minor  or  no  losses. 

In  the  QMF  approach,  a  signal  is  decomposed  into  two  frequency  bands,  namely 
below  and  above  f/4,  with  fs  being  the  sampling  frequency.  In  addition,  both  signal 
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components  are  downsampled  by  two,  thus  keeping  the  same  data  rate.  This  is  shown  in 
the  analysis  network  in  Figure  10.  In  the  synthesis  section  the  signal  is  reconstructed, 
ideally  without  losses  or  distortion.  Since  the  downsampling  operation  causes  aliasing, 
particular  care  has  to  be  taken  in  designing  the  two  pairs  of  lowpass  filter  (LPF)  and 
highpass  filter  (HPF). 

By  doing  so,  one  gets  two  half-sampled  signals  at  half  the  initial  bandwidth  each. 
Both  signals  are  now  downsampled  by  two  so  every  other  sample  is  dropped.  Finally,  for 
the  analysis  stage  one  gets  the  signal  yi(n),  which  corresponds  to  the  output  of  the 
lowpass  filter  and  the  signal  y2(n),  which  corresponds  to  the  highpass  filter  output.  Again, 
the  above  procedure  can  be  seen  in  Figure  10  below. 


Figure  10.  Two  Subband  QMF  Analysis-Synthesis  Scheme 


It  can  be  shown  that  one  of  the  conditions  for  perfect  reconstruction  is  that  the 
sum  of  the  Fourier  transforms  H L( e! 01  )  and  HH(eico)  of  hL(n)  and  hH(n),  respectively,  must 
be  at  every  frequency  equal  to  one.  [Jain] 


+  HH(eMy  =  1 


(1.13) 


Figure  1 1  shows  the  highpass  response  magnitude  as  the  mirror  image  of  the  lowpass 

7T 

response  magnitude  with  respect  to  the  middle  frequency  —  (Quadrature  frequency). 
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Figure  1 1 .  Quadrature  Frequencies 

One  of  the  characteristics  of  this  decomposition  is  that  it  preserves  the  energy  of 
the  signals,  in  the  sense  that 

EkiWI+Eh  ikf =ZkW|2  o-14) 

k  k  n 

with  yi,  y2  as  the  low  frequency  and  high  frequency  components.  Since  this 
decomposition  will  be  perfonned  recursively  in  all  components,  the  overall  energy  will 
be  preserved  throughout  the  overall  decomposition.  Likewise,  if  the  same  signal  x  is 
decomposed  into  N components  yi,  V2,  ■■■,  V/v,  the  above  fonnula  can  be  expanded  as 

ZZL>#]|2 =ZhH2  (U5) 

i=l  k  n 

The  reconstruction  process  is  very  similar  to  the  decomposition.  It  is  very  easy  to 
show  that  from  the  two  outputs  yi(n)  and  y2(n),  one  can  reconstruct  a  signal  x’(n)  in  a 
way  that  theoretically  is  identical  to  the  original  signal  x(n).  The  summary  of  the 
reconstruction  is  the  two  subband  signals  are  interpolated  by  a  factor  of  two  and  then 
after  filtering  them  with  another  pair  of  filters,  similarly  made  as  in  the  analysis,  the  final 
signal  x  ’( n)  is  produced.  [Swaminathan] 
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D.  DAUBECHIES  WAVELETS 


For  the  implementation  of  the  QMF  bank  a  class  of  low  complexity  digital  filters 
is  used,  namely  the  Daubechies  family,  that  satisfy  the  property  of  the  Quadrature  Mirror 
Filters  for  energy  preservation. 

With  the  Daubechies  filters  included  in  Matlab’s  “Wavelet  Toolbox,”  the  perfect 
reconstruction  condition  is  satisfied  with  FIR  filters  of  low  order.  The  command  for  these 
filters  is  “wfilters(‘dBN’),”  where  N  is  related  to  the  order.  Each  Daubechies  filter 
generates  a  corresponding  Daubechies  Wavelet,  shown  in  Figure  12.  In  this  family  of 
filters,  for  N=1  one  has  the  simplest  wavelet,  the  Haar  Wavelet. 


db2  db3  db4  db5  db6 


db7  db8  db9  dblO 


Figure  12.  Examples  of  the  Order  of  Daubechies  Filters  [From:  Matlab] 

If  Fo  is  a  Daubechies  scaling  filter  with  F )=QMF(F 'o)  for  its  QMF,  then  the 
transfer  functions  of  Fo  and  F]  will  be  according  to  Equation  1.12  [Matlab]: 

\F0(a>)\2 +\Fl(a>)\2  =1  (1.16) 

In  order  to  visualize  this,  one  can  follow  an  example  with  “dB10”  found  in  the 
Matlab’s  Wavelet  Toolbox.  On  the  first  two  plots  in  Figure  13,  there  is  the  lowpass  filter 
for  dB10  and  its  highpass  Quadrature  Mirror  Filter.  In  the  second  row,  the  FFTs  of  the 
above  filters,  respectively,  can  be  seen.  Finally,  in  the  last  row  there  is  the  sum  of  the 
absolute  values  of  the  squared  FFTs  that  are  equal  to  one  at  each  frequency,  and  thus 
proving  that  the  QMF  condition  is  satisfied  [Matlab]. 
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(1.17) 


|fft(filter)|~  +  fft(qmf(filter)) 2  =1 


dbl  0  low-pass  filter 
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Figure  13.  Example  of  the  QMF  Condition  for  the  dB10  [From:  Matlab] 


E.  THE  FLATTENING  OF  NOISE  IN  THE  SUBBANDS 

1.  Introduction 

Below,  it  will  be  seen  how  the  signal  decomposition  via  the  use  of  filter  hanks 
helps  to  better  estimate  the  frequencies  in  the  subbands  than  in  the  fullband  by  increasing 
the  local  SNR.  The  filters  used  are  approximately  ideal  and  there  is  minimal  overlap 
between  adjacent  filters.  It  is  well-known  that  parametric  techniques  to  estimate  signal 
models  do  not  work  well  in  the  presence  of  colored  noise.  This  is  particularly  true  when 
this  is  applied  to  spectral  estimation.  The  use  of  frequency  decomposition  partially 
circumvents  this  problem,  since  in  every  subband  the  noise  is  approximately  white.  When 
colored  noise  is  added  in  the  input  signal  with  the  individual  process  in  the  subbands,  it 
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can  be  assumed  to  be  white.  Nevertheless,  the  same  assumption  could  not  be  made  in  the 
fullband.  That  is  extremely  useful  for  this  case,  where  the  statistics  of  the  input  signal  not 
only  are  not  known  a  priori,  but  also  have  to  be  estimated  in  space  and  determined 
adaptively. 


2,  SNR  Reduction 


In  order  to  be  able  to  visualize  the  SNR  reduction  that  is  accomplished  with  the 
filter  bank  decomposition,  a  number  M  of  sinusoids  buried  in  complex  noise  n(n)  is 
assumed.  Then  the  signal  model  x(n)  becomes 

M- 1 

x(n)  =  ^  Aiej<a,n  +  n(n),  Ai  =  \  At  |  eJ (1-18) 

(=0 

where  the  amplitudes  A/(  are  unknown  but  constant  with  uniformly  distributed  phase 
angles  in  the  [-7t,jt).  For  the  noise,  it  is  assumed  to  be  a  zero  mean  wide  sense  stationary 
(WSS)  random  process  uncorrelated  with  the  sinusoids. 


Under  the  assumption  that  M  is  large  enough,  then  x(n)  is  a  zero  mean  WSS 
process  with  the  following  autocorrelation  function  [Tcacenko]. 


M- 1 

f;  =  | 4|  (1 19) 

i= 0 

If  the  local  signal-to-noise  ratio  is  defined  for  each  of  the  M  sinusoids  with  respect 
to  the  noise  as  SNRj,  then 


SNR  =  - 


i=0, 


To  keep  the  model  simple  ideal  filters  are  assumed.  Therefore, 

R^YMP.e-'+RJk),  I  A2™ 

CO:  G/ 


M 


M 


(1.20) 


(1.21) 


where Im  is  the  band  covered  by  the  mh  filter.  Combining  the  above  Equations  1.19,  1.20, 
1.21,  one  gets 


SNR,sub=M(SNRiJull )  (1.22) 

which  illustrates  that  for  each  of  the  i  sinusoids,  the  signal-to-noise  ratios  of  the  mth 
subband  is  greater  than  the  equivalent  one  in  the  whole  band  by  a  factor  of  M,  inversely 
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proportional  to  the  bandwidth  of  each  filter.  From  that,  it  is  easy  to  see  that  the  estimation 
in  the  subbands  is  going  to  be  more  accurate  than  in  the  fullband. 


At  this  point,  it  needs  to  be  mentioned  that  increasing  the  decimation  of  the  whole 
band  to  a  large  number  M,  in  practice,  is  not  always  as  profitable  as  one  might  think.  The 
reason  is  that  the  length  of  the  data  in  each  subband  will  be 


N,+Nf-2 

M 


+  1 


(1.23) 


where  Ns  is  the  number  of  available  observations  of  the  signal  x(n)  and  Nf  is  the  length  of 
each  of  the  analysis  filters.  From  here,  it  can  be  seen  that  the  brute  force  approach  of 
increasing  the  length  of  each  subband  would  result  in  the  increase  of  filter  length  and 
equitably  to  the  increase  of  complexity  of  the  filter  bank.  As  a  general  rule,  it  could  be 
said  that  the  most  appropriate  size  for  M  would  be  as  large  as  possible  so  there  would  be 
enough  samples  in  each  subband. 


3.  The  Whitening  of  Noise  in  the  Subbands 

It  is  assumed  that  the  noise  n(n)  of  the  signal  x(n)  has  a  nonconstant  power 
spectral  density  as  in  Figure  14,  and  the  magnitude  square  response  of  the  mtb  analysis 
filter  is  as  shown  in  Figure  15. 


Figure  14.  Nonconstant  Power  Spectral  Density  [From:  Tkacenko] 
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Figure  15.  Power  Spectrum  [From:  Tkacenko] 


It  is  very  easy  to  see  that  while  the  whole  noise  spectrum  [0,2n]  is  not  constant, 


the  smaller  area  inside  the  kth  subband 


Ink  2n(k  +  X) 
M’M 


is  approximately  constant.  In 


other  words,  the  larger  the  selected  M  is,  the  smaller  the  error  in  the  approximation  of 
colored  noise  with  white.  However,  it  should  be  kept  in  mind,  as  mentioned  above,  that 
M  can  not  be  arbitrarily  large. 


According  to  [Tkacenko],  with  the  usage  of  the  spectral  flatness  measure  y]  of  a 

WSS  random  process  and  based  on  three  properties  mentioned  below,  he  proves  that  the 
weighted  geometric  mean  of  the  spectral  flatness  values  in  the  subbands  is  greater  than  or 
equal  to  the  flatness  in  the  fullband.  The  three  properties  are  that  the  filter  bank  is 
maximally  decimated, 


M- 1  I 

I—i 

i=0  >h 


that  the  magnitude  squared  response  of  the  i,h  filter  is  Nyquist(/?/), 


Vz 


(1.24) 

(1.25) 


and  that  the  set  of  analyses  satisfies  a  generalized  version  of  the  power  complementary 


property, 
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=  1  Vo 


(1.26) 


n, 

The  review  of  the  whole  analysis  exceeds  the  objective  of  this  thesis  and  will  not 
be  covered  here.  [Tkacenko] 

F.  SUMMARY 

As  shown  in  the  background  theory  and  literature  review,  filter  hanks  can  be  very 
effective  in  the  estimation  of  frequencies  of  sinusoids  buried  in  noise.  In  particular,  when 
proper  filter  classes  are  used  for  the  implementation  of  the  FB,  the  SNR  in  the  subbands 
exceeds  that  in  the  fullband.  In  between  these  classes  are  the  QMF  and  the  FIR  filters 
which  the  author  chose  to  test  here.  For  both  of  them,  it  will  be  shown  how  the 
introduced  colored  noise  affects  their  behavior  and  therefore  the  measurements. 
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II.  SIGNAL  DECOMPOSITION  BY  FILTER  BANKS 


In  this  chapter  a  general  review  of  the  models  used  in  the  decomposition  of  the 
signal  using  filter  banks  is  presented.  In  the  first  approach,  Quadrature  Mirror  Filters 
filters  were  used  and  in  the  second  approach,  DFT  FB.  Each  case  was  implemented  in 
Matlab’s  Simulink  and  the  results  were  compared  in  terms  of  accurate  frequency 
estimation. 

A.  FLEXIBLE  STRUCTURE  SPACE  MODEL 

The  signal  used  is  generated  by  the  plant  presented  in  [Cristi],  modeling  the 
behavior  of  a  flexible  structure  space  model.  The  main  parts  are  the  rigid  body  and  the 
flexible  appendages.  The  rigid  body  is  modeled  as  a  double  integrator  ( I/s2)  and  the 
flexible  appendages  as  a  parallel  system  with  underdamped  poles.  A  proportional  integral 
derivative  (PID)  controller  is  also  implemented  in  a  closed  loop  to  track  an  external 
signal,  a  combination  of  a  step  function  and  a  sine  wave.  The  filter  bank  decomposition 
the  author  implemented  is  added  to  the  output  and  input  signals.  The  described  model  as 
used  by  the  author  in  Matlab’s  Simulink,  can  be  seen  below  in  Figure  16.  Also  in  Figure 
17  the  model’s  rotation  angle  in  degrees  can  be  seen. 
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Figure  16.  Flexible  Structure  Space  Model 
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Figure  17.  FSS  Rotation  Angle 


Figure  18  shows  the  output  (y)  of  the  model  and  Figure  19  illustrates  the  input  (u). 
In  both  cases,  the  peak  in  the  beginning  from  the  step  function  and  the  ongoing 
flexibilities  can  be  seen.  The  filter  bank  decomposition  designed  in  both  cases  of  the 
Quadrature  Mirror  Filters  and  the  DFT  FB  is  going  to  be  applied  only  to  the  output  (y), 
but  could  also  be  implemented  for  the  input  ( u )  as  well. 
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y  Output 


Figure  18.  “y”  Output  from  FFS  model 
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Figure  19.  “u”  Input  to  FSS  model 
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B.  ANALYSIS  WITH  QUADRATURE  MIRROR  FILTERS 

As  mentioned  above,  the  frequency  components  of  the  vibrations  were  estimated 
by  decomposing  the  observed  signals  into  a  number  of  channels.  In  this  section  the  author 
uses  QMF  decomposition. 

1.  Analysis  with  Quadrature  Mirror  Filters  without  Colored  Noise 

The  decomposition  is  based  on  the  following  shown  in  Figure  20. 


Figure  20.  Filter  Bank  for  Two  QMFs  in  Matlab 


The  signal  y  coming  from  the  workspace  is  the  output  (y)  from  the  FSS  model,  as 
a  sampled  sequence  y[n]=y(nTs)  with  Ts=l/Fs  as  the  sampling  interval.  The  LP(z)  and 
BPO(z)  filters  represent  the  lowpass  filter  ( LPF)  and  the  highpass  filter  ( HPF)  in  Figure 
21  respectively.  The  two  bandlimited  outputs  yi  and  y2  are  maximally  decimated  by  a 
factor  of  two  in  Matlab  to  give  the  two  corresponding  components  y i[m  / ]  and  y2[m2] 
shown  in  Figure  21  with  the  same  data  rate  as  the  original  signal. 


yfn] 


y,[m,] 

y2[m2] 


Figure  21.  General  Model  Filter  Bank  for  Two  QMFs  [After:  Cristi] 
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The  two  components  are  sampled  at  a  rate  F/2  and  the  lowpass  and  highpass 
components  of  the  frequency  spectrum  between  zero  and  F/2  are  mapped  as  shown  in 
Figure  22.  From  these  two  bands,  their  energy  is  calculated  by  taking  the  standard 
deviation  of  the  output.  For  the  band  that  has  the  bigger  energy,  the  choice  is  made  to 
decompose  it  again  as  before.  By  following  the  same  procedure,  the  signal  is  decomposed 
into  six  levels  and  seven  outputs  are  received  from  which  one  can  calculate  which  one 
has  finally  the  higher  energy.  The  decomposition  followed  for  the  FFS  model  output  is 
shown  in  Figure  23. 


4  2 


Figure  22.  Frequency  Bands  for  Two  QMFs  [After:  Cristi] 


Figure  23.  Frequency  Decomposition  Using  Cascaded  Filters 

Even  though  the  higher  energy  for  the  given  signal  is  always  going  to  be  in  the  LP 

band  because  of  the  step  function  used  in  the  FSS  model,  it  is  chosen  to  modify  the 

decomposition  after  the  third  layer  to  the  HP  band.  The  reason  is  to  estimate  the  dominant 
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frequency  by  the  higher  energy,  which  is  still  lower  than  the  energy  of  the  step  function. 
For  the  same  reason,  the  output  content  of  the  lower  frequency  component  caused  by  the 
step  function,  is  discarded. 

Based  on  the  Noble  Identity  shown  below  in  Figure24  [Cristi],  the  author  finally 
has  the  filter  bank  shown  in  Figure  25  with  the  frequency  decomposition  displayed  in 
Figure  26. 


Figure  25. 
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Figure  26.  Frequency  Bands  for  Eight  QMFs 

The  above  decomposition,  as  implemented  in  Matlab,  is  also  shown  in  Figure  27,  where 
if  g(z)  is  the  first  LP  filter  and  h(z)  is  the  first  HP  filter,  then  the  following  transfer 
functions  are  obtained:  LP(z)  =  g(z)g(z2)g(z4) ,  BPO(z)  =  g(z)g(z2)A(z4)g(z‘s), 
BP\(z)  =  g(z)g(z2)/2(z4)/?(z8)g(z16)g(z32) ,  BP2(z)  =  g(z)g(z2)h(z4)h(zs)g(zl6)h(zi2) , 
BP3(z)  =  g(z)g(z2)h(z4)h(z8)h(z16) ,  and  BP4(z)  =  g(z)h(z2)  and  LP(z )  =  h(z)  . 


Figure  27.  Filter  Bank  with  QMF  in  Matlab 
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For  Fs=lHz  the  following  two  pairs  of  graphs  (Figures  28  and  29)  with  tables 
(Tables  1  and  2)  are  obtained  that  show  the  energies  within  different  frequency  ranges. 
The  first  pair  shows  the  energy  of  the  outputs  for  dB15  and  the  second  one  for  dB45.  The 
rest  of  the  graphs  and  tables  for  two,  three,  four,  five,  six,  and  again  seven  outputs  for  the 
case  of  dB15  are  shown  in  Appendix  A.  At  this  point  it  needs  to  be  mentioned  that  the 
whole  analysis  the  author  followed  was  based  on  the  outcomes  of  the  dB15  case  and 
would  have  had  a  different  filter  bank  analysis  if  the  analysis  was  based  on  the  dB45 
case. 


128 


Frequency 


Figure  28.  Standard  Deviation  of  Seven  QMD  filters  for  dB  1 5 


Standard  Deviation  of  Seven  QMF  Filters  for  dB15 

Output 
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yi 

y2 

ys 

y4 
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y6 

Standard 

Deviation 

7.2592 

0.3434 

0.6996 

0.3848 

0.2636 

0.2747 

0.1019 

Table  1 .  Standard  Deviation  of  Seven  QMF  filters  for  dB  1 5 
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Figure  29.  Standard  Deviation  of  Seven  QMF  filters  for  dB45 


Standard  Deviation  of  Seven  QMF  filters  for  dB45 

Output 

yo 

yi 

y2 

ys 

y4 

ys 

ye 

Standard 

Deviation 

12.4913 

20.1186 

0.2717 

94.5267 

49.9890 

2.3487 

0.1494 

Table  2.  Standard  Deviation  of  Seven  QMF  filters  for  dB45 


Based  on  this  decomposition,  the  spectral  components  with  the  highest  energy  are 
known  and  the  frequency  response  of  the  system  within  this  frequency  band  can  be 
estimated.  Below  in  Figure  30  and  Table  3  the  wanted  frequency  subband  with  the  higher 
energy  is  isolated  with  the  use  of  weights. 
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Figure  30.  Standard  Deviation  of  Seven  Weighted  QMF  Filters  for  dB  1 5 


Standard  Deviation  of  Seven  weighted  QMF  filters  for  dB15 

Output 

yo 

yi 

y2 

ys 

y4 

ys 

yo 

Standard 

Deviation 

0 

0 

0.6996 

0.3848 

0 

0 

0 

Table  3.  Standard  Deviation  of  Seven  Weighted  QMF  Filters  for  dB15 

Furthermore,  in  order  for  the  system  to  be  adaptive  and  dynamic  the  author 
developed  a  system  of  loops  in  Matlab  that  selects  each  time  the  higher  output  of  the 
stage  and  decomposes  it  again.  The  loop  was  repeated  four  times  and  the  same  five 
outputs  as  in  the  case  from  the  previous  model  for  five  outputs  should  have  been 
obtained.  This  time  though,  the  selection  could  not  be  guided  like  was  the  case  before 
where  after  the  third  decomposition  the  author  chose  not  to  divide  again  the  higher-in- 
energy  subband  coming  from  the  step  response  in  the  lowpass  area.  Here  the  author  lets 
the  system  decide  for  every  decomposition  and  so  the  higher  energy  subband  is  as 
expected  in  the  lower  frequency  subband  with  different  values  for  the  graph  and  table 
than  before.  Bellow  in  Figure  31  is  the  model  used  for  that  dynamic  decomposition  and 
after  that  its  outputs  in  the  Tables  4  and  5  for  15  and  45dB,  respectively. 
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Figure  3 1 .  Dynamic  Filter  Bank  of  Four  stages 
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Standard  Deviation  of  Dynamic  Bank  with  Five  QMF  Filters  for  dB15 

Output 

yo 

yi 

y2 

ys 

y4 

Standard 

Deviation 

14.3204 

1.1609 

0.3977 

0.2747 

0.1019 

Table  4.  Standard  Deviation  of  Dynamic  Bank  with  Five  QMF  Filters  for  dB  15 


Standard  Deviation  of  Dynamic  Bank  with  Five  Filters  for  dB45 

Output 

yo 

yi 

y2 

y3 

y4 

Standard 

Deviation 

22.8628 

17.3112 

4.8914 

2.3487 

0.1494 

Table  5.  Standard  Deviation  of  Dynamic  Bank  with  Five  QMF  Filters  for  dB45 

Even  though  in  both  cases  of  dB15  and  dB45  the  higher-in-energy  output  can  be 
clearly  selected  each  time,  in  the  second  case  the  numbers  are  more  accurate  because  of 
the  higher  order  Daubechies  filter.  Finally  here,  weights  can  be  used  again  to  isolate  the 
outputs  that  are  preferred. 

2.  Analysis  with  Quadrature  Mirror  Filters  with  Colored  Noise 

For  the  author  to  be  able  to  fully  see  the  advantages  of  the  above  decomposition, 
colored  noise  is  added  to  the  model  as  shown  below  in  Figure  32,  with  the  power 
spectrum  shown  in  Figure  33. 
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Figure  32. 


FSS  Model  with  Colored  Noise  for  the  QMF 
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Figure  33.  Power  Spectrum  of  the  Used  Colored  Noise 

After  the  added  colored  noise,  the  standard  deviations  of  the  new  outputs  are 
again  calculated  and  the  author  gets  the  graphs  in  Figures  34  and  35  and  tables  in  Tables 
6  and  7.  From  these  values,  the  author  can  clearly  say  that  the  added  noise  did  not  affect 
the  model  since  the  values  are  analogous  to  the  previous  ones  when  the  model  had  no 
added  noise.  The  rest  of  the  graphs  and  tables  for  two,  three,  four,  five,  six  and  again 
seven  outputs  for  the  case  of  dB15  with  added  colored  noise  are  shown  in  Appendix  B. 


Power  Spectrum  of  colored  noise 
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Figure  34.  Power  Spectrum  with  Colored  Noise  of  Seven  QMF  Filters  for  dB  1 5 


Standard  Deviation  with  Colored  Noise  of  Seven  QMF  Filters  for  dB15 

Output 

yo 

yi 

y2 

ys 

y4 

ys 

ye 

Standard 

Deviation 

408370 

0.4769 

1.1954 

0.9899 

0.7239 

0.3727 

0.3545 

Table  6.  Standard  Deviation  with  Colored  Noise  of  Seven  QMF  Filters  for  dB15 


Standard  Deviation  with  Colored  Noise  of  Seven  QMF  Filters  for  dB45 

34 


Figure  35. 


Standard  Deviation  with  Colored  Noise  of  Seven  QMF  Filters  for  dB45 
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y3 
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Standard 

Deviation 

250980 

17.8080 

0.0321 

144.6512 

133.0226 

3.7517 

1.0448 

Table  7.  Standard  Deviation  with  Colored  Noise  of  Seven  QMF  Filters  for  dB45 


Analogous  results  are  gotten  if  the  same  colored  noise  is  applied  to  the  dynamic 
FB  for  five  outputs  (Tables  8  and  9). 


Standard  Deviation  with  Colored  Noise  of  Dynamic  Bank  with  Five  Filters  for  dB15 

Output 
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Deviation 

409440 

5.8021 

0.  7592 

0.  3727 

0.3545 

Table  8.  Standard  Deviation  with  Colored  Noise  of  Dynamic  Bank  with  Five  QMF  Filters 

for  dB15 


Standard  Deviation  with  Colored  Noise  of  Dynamic  Bank  with  Five  Filters  for  dB45 

Output 
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yi 

y2 
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y4 

Standard 

Deviation 

322720 

80.2683 

10.5877 

3.7517 

1.0448 

Table  9.  Standard  Deviation  with  Colored  Noise  of  Dynamic  Bank  with  Five  QMF  Filters 

for  dB45 

C.  ANALYSIS  WITH  DFT  FILTER  BANKS 

Another  class  of  filters  that  preserves  the  energy  of  a  signal  when  it  is  being 
decomposed  into  subbands  is  the  FIR  filters,  while  satisfying  the  conditions  described  in 
the  background  theory.  Even  though  the  author  was  not  able  to  use  ideal  filters  because  of 
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their  infinite  order,  an  approximation  was  made  via  the  use  of  a  FIR  filter  with  order  63 
as  the  original  filter.  Likewise,  the  order  can  very  easily  be  increased  for  more  accurate 
results,  which  goes  beyond  the  purpose  of  this  thesis.  Below  are  the  block  parameters  of 
the  original  filter  used  in  Matlab. 


Figure  36.  Block  Parameters  of  the  Used  FIR  Filter 


1.  Analysis  with  FIR  Filters  without  Colored  Noise 

For  the  analysis  with  the  above  filter,  a  filter  bank  was  created  using  the 
polyphase  components  in  Matlab  as  shown  in  Figure  37.  Each  one  of  the  new  filters  has 
64/8  coefficients  from  the  original  FIR  filter  and  processes  the  signal  with  a  delay  of zl 
from  the  previous  one.  From  that  filter  hank  comes  eight  outputs,  of  which  IFFT  as 
described  in  the  background  theory  is  taken.  Then,  the  standard  deviation  is  calculated  in 
each  channel.  From  the  values  obtained,  the  determination  is  made  as  to  which  subband 
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Figure  38.  Standard  Deviation  of  Eight  FIR  filters 


Standard  Deviation  of  Eight  FIR  filters 

Output 

Y0 
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y2 
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y4 

ys 

ye 

y? 

Standard 

Deviation 

0.0727 

0.0171 

0.0766 

0.0704 

0.0249 

0.0507 

0.0072 

0.037 

Table  10.  Standard  Deviation  of  Eight  FIR  filters 


2.  Analysis  with  FIR  Filters  with  Colored  Noise 

On  the  above  model  the  same  colored  noise  used  in  Section  III.B.2  for  the  QMF 
case  was  added.  The  results  for  the  colored  signal  can  be  seen  below  in  Figure  39  and 
Table  1 1.  If  these  results  are  compared  with  the  previous  ones,  namely  the  case  without 
the  colored  noise,  it  will  again  be  verified  that  they  are  analogous.  It  is  obvious  that  no 
matter  the  added  colored  noise,  the  filter  bank  was  able  to  distinguish  again,  with  high 
precision,  the  higher-in-energy  subbands. 
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Figure  39.  Standard  Deviation  of  Eight  FIR  Filters  with  Colored  Noise 


Standard  Deviation  of  Eight  FIR  Filters  with  Colored  Noise 
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11383.3 

1911.3 

12699.0 

11888.2 

1406.7 

7267.1 

6.4108 

5499.1 

Table  1 1 .  Standard  Deviation  of  Eight  FIR  Filters  with  Colored  Noise 


D.  SUMMARY 

In  this  chapter,  the  behavior  of  FB  implemented  with  QMF  and  FIR  filters  was 
examined.  For  each  case  the  author  applied  a  signal  initially  noise-free  and  then  with 
added  colored  noise.  The  measurements  taken  were  presented  in  tables  and  bar  charts  for 
easier  comparison.  All  the  findings  and  deductions  are  summarized  in  Chapter  III. 
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III.  CONCLUSIONS 


In  this  thesis  a  method  of  frequency  estimation  by  Quadrature  Mirror  Filters  and 
Filter  Banks  was  introduced.  They  both  provide  a  suitable  decomposition  of  the  signal  in 
tenns  of  the  various  frequency  components. 

Both  approaches  are  very  efficient  in  terms  of  the  number  of  operations.  The 
particular  advantage  is  that  they  provide  an  estimate  of  the  frequencies  within  a  specified 
interval. 

Indeed,  the  results  in  both  cases  satisfy  the  theoretical  expectations.  For  the  first 
case  with  the  QMF  it  was  seen  that  while  decomposing  the  original  signal  the  results  are 
leading  towards  the  same  subband.  Analogous  values  were  received  even  in  the  case  with 
the  dynamic  FB  for  five  outputs  with  the  mentioned  difference  due  to  the  altered 
decomposition  algorithm.  Likewise,  for  the  second  case  with  the  FIR  filters  the  standard 
deviation  in  all  the  subbands  increased  when  the  colored  noise  was  added,  but  still  in  an 
analogous  fashion. 

If  these  two  decomposition  Filter  Banks  are  compared,  it  can  be  seen  that  the 
QMF  for  the  same  number  of  outputs  gives  a  more  accurate  result,  meaning  that  the 
determined  subband  is  a  lot  narrower.  Moreover,  as  the  level  of  decomposition  is 
extended,  the  difference  between  them  also  grows. 

On  the  other  hand,  the  case  of  QMF  has  the  disadvantage  of  having  to  calculate 
the  energy  of  each  subband  on  each  level  of  the  decomposition.  That  procedure 
introduces  a  number  of  serious  calculations  which  results  in  an  increase  of  the  system’s 
feedback  time.  For  the  calculating  procedure  used  in  this  thesis,  the  FB  with  the  FIR 
filters  gave  the  faster  results  either  with  or  without  the  colored  noise. 

As  a  final  conclusion  it  could  be  said  that  both  cases  have  their  positive  and 
negative  sides,  so  the  selection  should  be  made  each  time  depending  on  the  environment 
and  requirements  of  the  current  application. 

Further  processing,  although  not  addressed  in  this  thesis,  is  to  refine  the  estimates 
by  properly  modeling  the  system  and  estimating  the  parameters.  In  this  case  the  two 
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proposed  approaches  have  the  advantage  that  they  use  the  signal  only  in  a  narrow  part  of 
the  spectrum,  thus  yielding  a  better  signal-to-noise  ratio  and  better  “whiteness”  of  the 
noise  spectrum.  Further  research  needs  to  be  undertaken  to  develop  this  method. 
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APPENDIX  A:  GRAPHS  AND  TABLES  OF  THE  QMF  FILTER 

BANK  FOR  DB15 


Figure  40.  Standard  Deviation  of  Two  QMF  Filters  for  dB15 


Standard  Deviation  of  Two  QMF  Filters  for  dB15 
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yi 
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1.5630 

0.1019 

Table  12.  Standard  Deviation  of  Two  QMF  Filters  for  dB15 


43 


10 


c 

o 

03 
> 
c D 
Q 

CO 


2: 


yo 


4 


Frequency 


7 

r 


2 


Figure  4 1 .  Standard  Deviation  of  Three  QMF  Filters  for  dB  1 5 


Standard  Deviation  of  Three  QMF  Filters  for  dB15 
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3.5030 

0.2747 

0.1019 

Table  13.  Standard  Deviation  of  Three  QMF  Filters  for  dB15 
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Figure  42.  Standard  Deviation  of  Four  QMF  Filters  for  dB  1 5 


Standard  Deviation  of  Four  QMF  Filters  for  dB15 
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7.2592 

0.3977 

0.2747 

0.1019 

Table  14.  Standard  Deviation  of  Four  QMF  Filters  for  dB15 
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Figure  43. 


Standard  Deviation  of  Five  QMF  Filters  for  dB15 


Standard  Deviation  of  Five  QMF  Filters  for  dB15 
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7.2592 

0.3434 

0.4345 

0.2747 

0.1019 

Table  15.  Standard  Deviation  of  Five  QMF  Filters  for  dB15 
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Figure  44.  Standard  Deviation  of  Six  QMF  Filters  for  dB  1 5 


Standard  Deviation  of  Six  Filters  for  dB15 
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7.2592 

0.3434 

0.5587 

0.2636 

0.2747 

0.1019 

Table  16.  Standard  Deviation  of  Six  QMF  Filters  for  dB15 
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Figure  45.  Standard  Deviation  of  Seven  QMF  Filters  for  dB15 


Standard  Deviation  of  Seven  QMF  Filters  for  dB15 
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7.2592 

0.3434 

0.6996 

0.3848 

0.2636 

0.2747 

0.1019 

Table  17.  Standard  Deviation  of  Seven  QMF  Filters  for  dB15 
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APPENDIX  B:  GRAPHS  AND  TABLES  OF  THE  QMF  FILTER 
BANK  FOR  DB15  WITH  ADDED  COLORED  NOISE 


Frequency 


Figure  46.  Standard  Deviation  of  Two  QMF  Filters  for  dB15  with  Added  Colored  Noise 


Standard  Deviation  of  Two  QMF  Filters  for  dB15  with  Added  Colored  Noise 
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274910 

0.3545 

Table  18.  Standard  Deviation  of  Two  QMF  filters  for  dB15  with  Added  Colored  Noise 
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Figure  47. 


Standard  Deviation  of  Three  QMF  Filters  for  dB15  with  Added  Colored  Noise 


Standard  Deviation  of  Three  QMF  Filters  for  dB15  with  Added  Colored  Noise 
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Table  19.  Standard  Deviation  of  Three  QMF  Filters  for  dB15  with  Added  Colored  Noise 
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Figure  48. 


Standard  Deviation  of  Four  QMF  filters  for  dB15  with  Added  Colored  Noise 


Standard  Deviation  of  Four  QMF  Filters  for  dB15  with  Added  Colored  Noise 
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408370 

0.  7592 

0.  3727 

0.3545 

Table  20.  Standard  Deviation  of  Four  QMF  Filters  for  dB15  with  Added  Colored  Noise 
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Figure  49. 


Standard  Deviation  of  Five  QMF  Filters  for  dB15  with  Added  Colored  Noise 


Standard  Deviation  of  Five  QMF  filters  for  dB15  with  Added  Colored  Noise 
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408370 

0.  4769 

0. 9357 

0.  3727 

0.3545 

Table  2 1 .  Standard  Deviation  of  Five  QMF  Filters  for  dB  1 5  with  Added  Colored  Noise 
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Figure  50.  Standard  Deviation  of  Six  QMF  Filters  for  dB  1 5  with  Added  Colored  Noise 


Standard  Deviation  of  six  Filters  for  dB15  with  Added  Colored  Noise 
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Table  22.  Standard  Deviation  of  Six  QMF  Filters  for  dB15  with  Added  Colored  Noise 
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Figure  5 1 .  Standard  Deviation  of  Seven  QMF  Filters  for  dB  1 5  with  Added  Colored 

Noise 


Standard  Deviation  of  Seven  QMF  Filters  for  dB15  with  Added  Colored  Noise 
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0.4769 

1.1954 

0.9899 

0.7239 

0.3727 

0.3545 

Table  23.  Standard  Deviation  of  Seven  QMF  Filters  for  dB15  with  Added  Colored  Noise 


54 


LIST  OF  REFERENCES 


[Cristi] 

[Cristil] 

[Galand] 

[Jain] 

[Karapatakis] 

[Lin] 

[Losada] 

[Manos] 

[Sandrock] 

[Swaminathan] 

[Tkacenko] 


R.  Cristi,  Control  of  flexible  Structures:  Frequency  Estimation, 
Unpublished  Report,  Department  of  Electrical  and  Computer 
Engineering,  Naval  Postgraduate  School,  Monterey,  California, 
September  2006. 

R.  Cristi,  Modern  Digital  Signal  Processing,  1st  edition,  Thomson 
Brooks/Cole,  Pacific  Grove,  California,  2004. 

Claude  R.  Galand  and  Henri  J.  Nussbaumer  “New  Quadrature  Mirror 
Filter  Structures,”  IEEE  Transactions  on  Acoustics,  Speech,  and 
Signal  Processing,  Vol.  ASSP-32,  No.  3,  June  1984. 

Vijay  K.  Jain  and  Ronald  E.  Crochiere,  “Quadrature  Mirror  Filter 
Design  in  the  Time  Domain,”  IEEE  Transactions  on  Acoustics, 
Speech,  and  Signal  Processing,  Vol.  ASSP-32,  No.  2,  April  1984. 

K.  Karapatakis,  Adaptive  Control  Of  A  Laser  Tracking  Device  For 
Space  Based  Application,  Master’s  Thesis,  Naval  Postgraduate 
School,,  Monterey,  California,  September  2006. 

Yuan-Pei  Lin  and  P.P.  Vaidyanathan,  “A  Kaiser  Window  Approach 
for  the  Design  of  Prototype  Filters  of  Cosine  Modulated  Filter  Banks,” 
IEEE  Signal  Processing  Letters,  Vol.  5,  No.  6,  June  1998. 

Ricardo  A.  Losada,  “Practical  FIR  Filter  Design  in  MATLAB,” 
Revision  1.1,  The  MathWorks,  Inc.,  3Apple  Hill  Dr.  Natick,  MA 
01760,  USA,  January  2004. 

E.  Manos,  Control  of  a  System  in  the  Presence  of  Flexible  Modes, 
Master’s  Thesis,  Expected  Completion  June  2007. 

Malte  Sandrock  and  Stefan  Schmitt,  “Realization  of  an  Adaptive 
Algorithm  with  Subband  Filtering  Approach  for  Acoustic  Echo 
Cancellation  in  Telecommunication  Applications,”  DSPesialists 
GmbH,  D- 10245,  Gennany. 

Kumar  Swaminathan  and  P.  P.  Vaidyanathan,  “Theory  and  Design  of 
Uniform  DFT,  Parallel,  Quadrature  Mirror  Filter  Banks,”  IEEE 
Transactions  on  Circuits  and  Systems,  Vol.  CAS-33,  No.  12, 
December  1986. 

Andre  Tkacenko  and  P.P.  Vaidyanathan,  “The  Role  of  Filter  Banks  in 
sinusoidal  Frequency  Estimation,”  California  Institute  of  Technology, 
136-93  Moore,  Pasadena,  California  91125. 


55 


[Vaidyanathan] 


Andre  Tkacenko  and  P.P.  Vaidyanathan,  “Sinusoidal  Frequency 
Estimation  Using  Filter  Banks,”  IEEE  Proceedings,  p.p.  3089-3092, 
2001. 


56 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Ft.  Belvoir,  Virginia 

2.  Dudley  Knox  Library 
Naval  Postgraduate  School 
Monterey,  California 

3.  Chairman,  Code  EC 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  California 

4.  Roberto  Cristi 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  California 

5.  Xiaoping  Yun 

Department  of  Electrical  and  Computer  Engineering 
Naval  Postgraduate  School 
Monterey,  California 

6.  Konstantinos  Tzellos 
Hellenic  Navy 
Athens,  Greece 


57 


